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Abstract The purpose of this brief comunication is to improve the hypothesis of the previous work 
of the author (de Divitiis, Theor Comput Fluid Dyn, doi:10.1007/s00162-010-0211-9) dealing with the 
finite-scale Lyapunov analysis of isotropic turbulence. There, the analytical expression of the structure 
function of the longitudinal velocity difference Au r is derived through a statistical analysis of the 
Fourier transformed Navier-Stokes equations, and by means of considerations regarding the scales of 
the velocity fluctuations, which arise from the Kolmogorov theory. Due to these latter considerations, 

^J ' this Lyapunov analysis seems to need some of the results of the Kolmogorov theory. 

This work proposes a more rigorous demonstration which leads to the same structure function, 
without using the Kolmogorov scale. This proof assumes that pair and triple longitudinal correlations 

^j , are sufficient to determine the statistics of Au r , and adopts a reasonable canonical decomposition 

of the velocity difference in terms of proper stochastic variables which are adequate to describe the 

^L | mechanism of kinetic energy cascade. 
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1 Introduction 



> 

o 

In the previous work [2], the author applies the finite-scale Lyapunov theory to analyse the homoge- 
l/~j , neous isotropic turbulence. In particular, this theory leads to the analytical closure of the von Karman- 

-^J • Howarth equation, giving the longitudinal triple velocity correlation k in terms of the longitudinal 

velocity correlation / and df/dr (see also the Appendix), and shows that the structure function of the 
£\« . longitudinal velocity difference is expressed by 

Au r = e + ^(x(^ 2 -i)-(c 2 -i)) (l) 

VIP^F) Vl + 2V> 2 (1 + X 2 ) 

where Au r = (u(x + r) — u(x)) ■ r/r, f , r\ and £ are uncorrelated centered gaussian random variables, 
with (£ 2 ) = (rj 1 ) = (C 2 ) = 1, and ip is a function of the Taylor scale Reynolds number R\ = u\t/v 
and of the separation distance r = |r| 



*"•*»- « M (2) 
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being u — \/ {v%.) the longitudinal velocity standard deviation, At is the Taylor microscale, and \ = 
x(-Ra) 7^ 1 is a proper function of R\ which provides nonzero skewness of Au r [2], [3]. 

In Ref. [2], the demonstration of Eq. ([TJ is carried out through statistical elements regarding the 
Fourier transformed Navier-Stokes equations, whereas the proof of Eq. @ is based on the fact that, 
according to the Kolmogorov theory, the ratio (small scale velocity)-(large scale velocity) depends on 
Xt/£ ~ \/R\, where £ is the Kolmogorov microscale. Therefore, the analysis of Ref. [2J seems to require 
the adoption of £ whose definition is based on another theory. 

Here, instead of using the Kolmogorov scale, we obtain Eqs. ([1]) and (J2J), starting from the canonical 
decomposition of the fluid velocity fluctuation in terms of proper centered random variables £& . In order 
to describe the mechanism of kinetic energy cascade, the variables £& are properly chosen in such a way 
that each of them exhibits a non-symmetric distribution function, with |(£f)| >>> (£fc) 3 ^ 2 : (£|) 3/ ' 4 
(k = 1, 2, ...). Moreover, we postulate that the knowledge of / and k represents a sufficient condition 
to determine the statistics of Au r . 



2 Lyapunov analysis of the velocity fluctuations 

This section renews the procedure for calculating the velocity fluctuations, which is based on the 
Lyapunov analysis of the fluid strain [2], and on the momentum Navier-Stokes equations 

du k du k , 1 dT kh 

~^r = ~^r~ Uh ^ — « — ( 3 

at axh p oxh 

being u = (ui, U2, 1*3), Tkh and p the fluid velocity, stress tensor and density, respectively. 

In order to obtain the velocity fluctuation, consider now the relative motion between two contiguous 
particles, expressed by the infinitesimal separation vector dx which obeys to the equation 

dx = Vu dx (4) 

where dx varies according to the velocity gradient which in turn follows the Navier-Stokes equations. 
As observed in Ref. [2], dx is much faster than the fluid state variables, and the Lyapunov analysis 
of Eq. (|2|) provides the expression of the local deformation in terms of maximal Lyapunov exponent 
A > 

dx 
dx 



e A(t-to) ( 5 ) 



The map x : Xo -> x, is the function which gives the current position x of a fluid particle located at 
the referential position xo at t = to 8 . Equation (J3j) can be written in terms of the referential position 

x m 

du k ( du k , ldT kh \ dx Qp 



dt \ 8xq p p 8xq p J dxh 

The adoption of the referential coordinates allows to factorize the velocity fluctuation and to express 
it in Lyapunov exponential form of the local deformation. As this deformation is assumed to be much 
more rapid than —duk/dxo P Uh + l/p dTkh/dxo p , the velocity fluctuation can be obtained integrating 
Eq. ([6]) with respect to the time, where —duk/dxo P Uh + 1/p dTkh / 9xq p is considered to be constant 

1 / du k , ldT kh \ 
A \ dx p p ox p J t=to 

This assumption is justified by the fact that, according to the classical formulation of motion of 
continuum media [8], the terms into the circular brackets of Eq. ([6]) are considered to be smooth 
functions of t -at least during the period of a fluctuation- whereas the fluid deformation varies very 
rapidly according to Eqs. ((5)-©. 
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3 Statistical analysis of velocity difference 

As explained in this section, the Lyapunov analysis of the local deformation and some plausible as- 
sumptions about the statistics of u lead to determine the structure function of Au r and its PDF. 

The statistical properties of Au r are here investigated expressing the velocity fluctuation through 
the following canonical decomposition [5] 

u = ^U fe 6, (8) 

k 

where Ufc {k = 1,2, ...) are proper coordinate functions of t and x, and £& (k — 1,2, ...) are certain 
dimensionless independent stochastic variables which satisfy 

<a) = o, (te j ) = 6 ij , (tejtk)=m jk P, W>»1, (#) = 0(i) (9) 

where Wijk = 1 for i = j = k, else tUijk=0- It is worth to remark that the variables £& are adequately 
chosen in such a way that they describe properly the mechanism of energy cascade. This latter property 
is expressed by the condition |(£j?)| = \p\ >>> 1. 

The dimensionless velocity fluctuation u is obtained in terms of £& substituting Eq. (JSJ) into Eq. ((7]) 

Uh = ^2 Ahij£i€j + -n-^] b h kik (10) 

ij k 

where r = f At and Uh = ii-hU- Therefore ^ li ■ Ahij£,i£,j arises from the inertia and pressure terms, 
whereas 1/ R\ ^2 k bhk^k is due to the fluid viscosity. Now, thanks to the local isotropy, Uh is a gaussian 
stochastic variable [5] , [5] , accordingly, £& satisfy into Eq. (TTU|) , the Lindeberg condition, a very general 
necessary and sufficient condition for satisfying the central limit theorem [5|. This condition does 
not apply to the velocity difference. In fact, as An is the difference between two correlated gaussian 
variables, its PDF could be a non gaussian distribution function. To study this question, the fluctuation 
Au r is first expressed in terms of £& 

Au r {r) = Y, AA riJ (r)^ + i- ^ Ab rk (r)£ k = L + S + G + + G~ (11) 

ij k 

This fluctuation can be reduced to the contributions L, S, G + and G~, appearing into Eq. (fTI) [B]: 
in particular, L is the sum of all linear terms due to the fluid viscosity S = /Sy&£j is the sum of 
all bilinear forms arising from the inertia and pressure terms, whereas G + and G~ are, respectively, 
definite positive and negative quadratic forms of centered gaussian variables, which derive from the 
inertia and pressure terms. The quantity L + S tends to a gaussian random variable being the sum of 
statistically orthogonal terms [B] , [5] , while G + and G~ are determined by means of the hypotheses of 
isotropy and of fully developed flow 

G-=-(( 2 -l)Mr), 

(12) 
G+ = (rj 2 - l)^(r) 

Observe that, due to these hypotheses, G + and G~ are uncorrected, thus 77, ( are two independent 
centered gaussian variables, with {rj 2 ) = (C 2 ) =1- Furthermore, as the knowledge of / and k is consid- 
ered to be a sufficient condition for determining the statistics of Au r , ip2 and i/^ are proportional with 
each other through a constant which depends only on R\ 

x{Rx) = v^w (13) 

Therefore, the longitudinal velocity difference can be written as 

Au r = i^ + i, 2 ( X ( V 2 - 1) - (C 2 - 1)) (14) 

where £ is a centered gaussian random variable with (£ 2 ) — 1, that thanks to the hypotheses of fully 
developed flow and of isotropy, is considered to be statistically independent from 77 and £. 
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Comparing the terms of Eqs. (J14I) and (TTTj) . we obtain that i/'i and ip2 are related with each other 
and that their ratio ip = ipi/ip2 depends on R\ and r 

?M ((G+ + G-) 2 ) , s 

Now, the divisor at the R. H. S. of Eq. (|15p is the sum of the following three terms: 

A = S ij S P9 (t i Z J t P Z q ), B = jLsijite^AWk, C = ^AblM) (16) 

Rx R x 

Hence, taking into account the properties © of £ fc (|(£f)| >>> 1, (££) = 0(1)), \B\ >» \A\, \C\, 
thus ip tends to a quantity arising only from the terms (£|) which appear in Eq. (|15[) . 

^ = tl = lp (r) x /R~ x (17) 

This expression corresponds to that obtained in Ref. [2] 



Mr, Rx) = J^fE ^ (r) ' ^ (0) = ° (1) (18) 

and the dimensionless longitudinal velocity difference is given by Eq. (JTJ) 
Au r _e + ^(x(^ 2 "l)-(C 2 -l)) 



J^A^f) v/l + 2^ 2 (l + X 2 ) 



(19) 



It is worth to remark that ip expresses the fluctuations ratio (large scale velocity)-(small scale velocity), 
that is ip ps u/u s ps (u 2 /\t){I s /u 2 ), where Z s and u s are, respectively the characteristic small scale 
and the corresponding velocity. This implies that u/u s ~ Xt/I s ~ y/R\, thus Z s is of the order of the 
Kolmogorov scale and u s l s /v « 1 is the corresponding Reynolds number. 

The distribution function of Au r is then expressed through the Frobenius-Perron equation [7], 
taking into account that £, r\ and £ are independent stochastic variables 

F(Au' r ) = [ [ [p(Op(v)p(0 S(Au' r -Au r (^ V ,0)d^ dr) d( (20) 

where Au r (£, i], £) is determined with Eq. (fT5|) . <5 is the Dirac delta and p is a centered gaussian PDF 
with standard deviation equal to the unity. The dimensionless statistical moments of Au r are easily 
calculated 

„ ((Au r ) n ) 



—2 E (l) ^ k (C~ k )((x(v 2 i) - (C 2 - i)) fe > (21) 



((Au r )T /2 (l + 2(l + x 2 )^ 2 ) , (| 
In particular, if 3, related to the mechanism of energy cascade, is 

8V> 3 (x 3 ~ 1) ,__* 

31 >~ (1 + 2^(1 +x 2))3/2 ^ 

In conclusion, x = x(R\) is implicitly calculated, in function of "0(0), taking into account that this 
Lyapunov theory gives #3(0) = —3/7 (see Appendix) [5] 

wm - 8^ 3 (0)(X 3 -1) _ 3 

31 j " (1 + 2V> 2 (0)(1 + X 2 ))3/ 2 " 7 [ ' 



where ip(0) — 1.075 is estimated as in Ref. [2J, and -(/((ZW,.) 2 ) and V>(r) are calculated in function of 
/(r) and fc(r). 
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5 Appendix 

For sake of convenience, this section reports some of the results dealing with the closure of the von 
Karman-Howarth equation, obtained in Refs. [5] and [3]. 

For fully developed isotropic homogeneous turbulence, the pair correlation function 
= K(xK(x + r)) 
u 2 
satisfies the von Karman-Howarth equation [4] 

df K(r) fd 2 f 4 ©A 101/ 

the boundary conditions of which are 

/(0) = 1, 

(26) 
lim /(r) = 



where u = -\/(?i 2 (x)) satishes the equation of the turbulent kinetic energy [3] 

-*- = - aT (27) 



and At = y/—l/f"(0) is the Taylor scale. This equation, obtained putting r = into Eq. (|2"5)) . gives 
the rate of kinetic energy in function of u and At [I], OP- The function K{r), giving the mechanism of 
energy cascade, is related to the longitudinal triple velocity correlation function k 



.3/ 9 , 4\,^ „ u __ lu _^(xK(x + r)) 



K(r) = u \w + - fc M, where fc ( r ) = o — ^^^ (28) 

\or r J u 6 

Thus, the von Karman-Howarth equation provides the relationship between the statistical moments 
((Z\w r ) 2 ) and ((ziu r ) 3 ) in function of r. 

The Lyapunov theory proposed in Ref. [5] gives the closure of the von Karman-Howarth equation, 
and expresses K(r) in terms of / and df/dr 



3 /!-/ df 



^-VY^ (29) 

Accordingly, the skewness oiAu r is [1] 

((/W) 3 ) _ 6fc(r) 

" 3lrJ " ((/A Ur )2)3/2 (2(1 - /(r))) 3 / 2 ldUj 

Therefore, the skewness of du r /dr is 

#3(0) = — (31) 
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